A Hybrid-parallel Architecture for Applications in Bioinformatics

Since the advent of Next Generation Sequencing (NGS) technology, the amount of data from whole genome sequencing has been rising fast. In turn, the availability of these resources led to the tapping of whole new research fields in molecular and cellular biology, producing even more data. On the other hand, the available computational power is only increasing linearly. In recent years though, special-purpose high-performance devices started to become prevalent in today’s scientific data centers, namely graphics processing units (GPUs) and, to a lesser extent, field-programmable gate arrays (FPGAs). Driven by the need for performance, developers started porting regular applications to GPU frameworks and FPGA configurations to exploit the special operations only these devices may perform in a timely manner. However, applications using both accelerator technologies are still rare. Major challenges in joint GPU/FPGA application development include the required deep knowledge of associated programming paradigms and the efficient communication both types of devices. In this work, two algorithms from bioinformatics are implemented on a custom hybrid-parallel hardware architecture and a highly concurrent software platform. It is shown that such a solution is not only possible to develop but also its ability to outperform implementations on similar- sized GPU or FPGA clusters in terms of both performance and energy consumption. Both algorithms analyze case/control data from genome- wide association studies to find interactions between two or three genes with different methods. Especially in the latter case, the newly available calculation power and method enables analyses of large data sets for the first time without occupying whole data centers for weeks. The success of the hybrid-parallel architecture proposal led to the development of a high- end array of FPGA/GPU accelerator pairs to provide even better runtimes and more possibilities

Large-scale genome-wide association studies on a GPU cluster using a CUDA-accelerated PGAS programming model

[Abstract] Detecting epistasis, such as 2-SNP interactions, in genome-wide association studies (GWAS) is an important but time consuming operation. Consequently, GPUs have already been used to accelerate these studies, reducing the runtime for moderately-sized datasets to less than 1 hour. However, single-GPU approaches cannot perform large-scale GWAS in reasonable time. In this work we present multiEpistSearch, a tool to detect epistasis that works on GPU clusters. While CUDA is used for parallelization within each GPU, the workload distribution among GPUs is performed with Unified Parallel C++ (UPC++), a novel extension of C++ that follows the Partitioned Global Address Space (PGAS) model. multiEpistSearch is able to analyze large-scale datasets with 5 million SNPs from 10,000 individuals in less than 3 hours using 24 NVIDIA GTX Titans.London. Wellcome Trust; 076113London. Wellcome Trust; 08547

FPGA-based Acceleration of Detecting Statistical Epistasis in GWAS

AbstractGenotype-by-genotype interactions (epistasis) are believed to be a significant source of unexplained genetic variation causing complex chronic diseases but have been ignored in genome-wide association studies (GWAS) due to the computational burden of analysis. In this work we show how to benefit from FPGA technology for highly parallel creation of contingency tables in a systolic chain with a subsequent statistical test. We present the implementation for the FPGA-based hardware platform RIVYERA S6-LX150 containing 128 Xilinx Spartan6-LX150 FPGAs. For performance evaluation we compare against the method iLOCi[9]. iLOCi claims to outperform other available tools in terms of accuracy. However, analysis of a dataset from the Wellcome Trust Case Control Consortium (WTCCC) with about 500,000 SNPs and 5,000 samples still takes about 19hours on a MacPro workstation with two Intel Xeon quad-core CPUs, while our FPGA-based implementation requires only 4minutes

Parallelizing Epistasis Detection in GWAS on FPGA and GPU-Accelerated Computing Systems

This is a post-peer-review, pre-copyedit version of an article published in IEEE - ACM Transactions on Computational Biology and Bioinformatics. The final authenticated version is available online at: http://dx.doi.org/10.1109/TCBB.2015.2389958[Abstract] High-throughput genotyping technologies (such as SNP-arrays) allow the rapid collection of up to a few million genetic markers of an individual. Detecting epistasis (based on 2-SNP interactions) in Genome-Wide Association Studies is an important but time consuming operation since statistical computations have to be performed for each pair of measured markers. Computational methods to detect epistasis therefore suffer from prohibitively long runtimes; e.g., processing a moderately-sized dataset consisting of about 500,000 SNPs and 5,000 samples requires several days using state-of-the-art tools on a standard 3 GHz CPU. In this paper, we demonstrate how this task can be accelerated using a combination of fine-grained and coarse-grained parallelism on two different computing systems. The first architecture is based on reconfigurable hardware (FPGAs) while the second architecture uses multiple GPUs connected to the same host. We show that both systems can achieve speedups of around four orders-of-magnitude compared to the sequential implementation. This significantly reduces the runtimes for detecting epistasis to only a few minutes for moderatelysized datasets and to a few hours for large-scale datasets.London. Wellcome Trust; 076113London. Wellcome Trust; 08547

Detailed stratified GWAS analysis for severe COVID-19 in four European populations

Given the highly variable clinical phenotype of Coronavirus disease 2019 (COVID-19), a deeper analysis of the host genetic contribution to severe COVID-19 is important to improve our understanding of underlying disease mechanisms. Here, we describe an extended genome-wide association meta-analysis of a well-characterized cohort of 3255 COVID-19 patients with respiratory failure and 12 488 population controls from Italy, Spain, Norway and Germany/Austria, including stratified analyses based on age, sex and disease severity, as well as targeted analyses of chromosome Y haplotypes, the human leukocyte antigen region and the SARS-CoV-2 peptidome. By inversion imputation, we traced a reported association at 17q21.31 to a ~0.9-Mb inversion polymorphism that creates two highly differentiated haplotypes and characterized the potential effects of the inversion in detail. Our data, together with the 5th release of summary statistics from the COVID-19 Host Genetics Initiative including non-Caucasian individuals, also identified a new locus at 19q13.33, including NAPSA, a gene which is expressed primarily in alveolar cells responsible for gas exchange in the lung.S.E.H. and C.A.S. partially supported genotyping through a philanthropic donation. A.F. and D.E. were supported by a grant from the German Federal Ministry of Education and COVID-19 grant Research (BMBF; ID:01KI20197); A.F., D.E. and F.D. were supported by the Deutsche Forschungsgemeinschaft Cluster of Excellence ‘Precision Medicine in Chronic Inflammation’ (EXC2167). D.E. was supported by the German Federal Ministry of Education and Research (BMBF) within the framework of the Computational Life Sciences funding concept (CompLS grant 031L0165). D.E., K.B. and S.B. acknowledge the Novo Nordisk Foundation (NNF14CC0001 and NNF17OC0027594). T.L.L., A.T. and O.Ö. were funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation), project numbers 279645989; 433116033; 437857095. M.W. and H.E. are supported by the German Research Foundation (DFG) through the Research Training Group 1743, ‘Genes, Environment and Inflammation’. L.V. received funding from: Ricerca Finalizzata Ministero della Salute (RF-2016-02364358), Italian Ministry of Health ‘CV PREVITAL’—strategie di prevenzione primaria cardiovascolare primaria nella popolazione italiana; The European Union (EU) Programme Horizon 2020 (under grant agreement No. 777377) for the project LITMUS- and for the project ‘REVEAL’; Fondazione IRCCS Ca’ Granda ‘Ricerca corrente’, Fondazione Sviluppo Ca’ Granda ‘Liver-BIBLE’ (PR-0391), Fondazione IRCCS Ca’ Granda ‘5permille’ ‘COVID-19 Biobank’ (RC100017A). A.B. was supported by a grant from Fondazione Cariplo to Fondazione Tettamanti: ‘Bio-banking of Covid-19 patient samples to support national and international research (Covid-Bank). This research was partly funded by an MIUR grant to the Department of Medical Sciences, under the program ‘Dipartimenti di Eccellenza 2018–2022’. This study makes use of data generated by the GCAT-Genomes for Life. Cohort study of the Genomes of Catalonia, Fundació IGTP (The Institute for Health Science Research Germans Trias i Pujol) IGTP is part of the CERCA Program/Generalitat de Catalunya. GCAT is supported by Acción de Dinamización del ISCIII-MINECO and the Ministry of Health of the Generalitat of Catalunya (ADE 10/00026); the Agència de Gestió d’Ajuts Universitaris i de Recerca (AGAUR) (2017-SGR 529). M.M. received research funding from grant PI19/00335 Acción Estratégica en Salud, integrated in the Spanish National RDI Plan and financed by ISCIII-Subdirección General de Evaluación and the Fondo Europeo de Desarrollo Regional (European Regional Development Fund (FEDER)-Una manera de hacer Europa’). B.C. is supported by national grants PI18/01512. X.F. is supported by the VEIS project (001-P-001647) (co-funded by the European Regional Development Fund (ERDF), ‘A way to build Europe’). Additional data included in this study were obtained in part by the COVICAT Study Group (Cohort Covid de Catalunya) supported by IsGlobal and IGTP, European Institute of Innovation & Technology (EIT), a body of the European Union, COVID-19 Rapid Response activity 73A and SR20-01024 La Caixa Foundation. A.J. and S.M. were supported by the Spanish Ministry of Economy and Competitiveness (grant numbers: PSE-010000-2006-6 and IPT-010000-2010-36). A.J. was also supported by national grant PI17/00019 from the Acción Estratégica en Salud (ISCIII) and the European Regional Development Fund (FEDER). The Basque Biobank, a hospital-related platform that also involves all Osakidetza health centres, the Basque government’s Department of Health and Onkologikoa, is operated by the Basque Foundation for Health Innovation and Research-BIOEF. M.C. received Grants BFU2016-77244-R and PID2019-107836RB-I00 funded by the Agencia Estatal de Investigación (AEI, Spain) and the European Regional Development Fund (FEDER, EU). M.R.G., J.A.H., R.G.D. and D.M.M. are supported by the ‘Spanish Ministry of Economy, Innovation and Competition, the Instituto de Salud Carlos III’ (PI19/01404, PI16/01842, PI19/00589, PI17/00535 and GLD19/00100) and by the Andalussian government (Proyectos Estratégicos-Fondos Feder PE-0451-2018, COVID-Premed, COVID GWAs). The position held by Itziar de Rojas Salarich is funded by grant FI20/00215, PFIS Contratos Predoctorales de Formación en Investigación en Salud. Enrique Calderón’s team is supported by CIBER of Epidemiology and Public Health (CIBERESP), ‘Instituto de Salud Carlos III’. J.C.H. reports grants from Research Council of Norway grant no 312780 during the conduct of the study. E.S. reports grants from Research Council of Norway grant no. 312769. The BioMaterialBank Nord is supported by the German Center for Lung Research (DZL), Airway Research Center North (ARCN). The BioMaterialBank Nord is member of popgen 2.0 network (P2N). P.K. Bergisch Gladbach, Germany and the Cologne Excellence Cluster on Cellular Stress Responses in Aging-Associated Diseases, University of Cologne, Cologne, Germany. He is supported by the German Federal Ministry of Education and Research (BMBF). O.A.C. is supported by the German Federal Ministry of Research and Education and is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy—CECAD, EXC 2030–390661388. The COMRI cohort is funded by Technical University of Munich, Munich, Germany. This work was supported by grants of the Rolf M. Schwiete Stiftung, the Saarland University, BMBF and The States of Saarland and Lower Saxony. K.U.L. is supported by the German Research Foundation (DFG, LU-1944/3-1). Genotyping for the BoSCO study is funded by the Institute of Human Genetics, University Hospital Bonn. F.H. was supported by the Bavarian State Ministry for Science and Arts. Part of the genotyping was supported by a grant to A.R. from the German Federal Ministry of Education and Research (BMBF, grant: 01ED1619A, European Alzheimer DNA BioBank, EADB) within the context of the EU Joint Programme—Neurodegenerative Disease Research (JPND). Additional funding was derived from the German Research Foundation (DFG) grant: RA 1971/6-1 to A.R. P.R. is supported by the DFG (CCGA Sequencing Centre and DFG ExC2167 PMI and by SH state funds for COVID19 research). F.T. is supported by the Clinician Scientist Program of the Deutsche Forschungsgemeinschaft Cluster of Excellence ‘Precision Medicine in Chronic Inflammation’ (EXC2167). C.L. and J.H. are supported by the German Center for Infection Research (DZIF). T.B., M.M.B., O.W. und A.H. are supported by the Stiftung Universitätsmedizin Essen. M.A.-H. was supported by Juan de la Cierva Incorporacion program, grant IJC2018-035131-I funded by MCIN/AEI/10.13039/501100011033. E.C.S. is supported by the Deutsche Forschungsgemeinschaft (DFG; SCHU 2419/2-1).Peer reviewe

Detailed stratified GWAS analysis for severe COVID-19 in four European populations

Given the highly variable clinical phenotype of Coronavirus disease 2019 (COVID-19), a deeper analysis of the host genetic contribution to severe COVID-19 is important to improve our understanding of underlying disease mechanisms. Here, we describe an extended GWAS meta-analysis of a well-characterized cohort of 3,260 COVID-19 patients with respiratory failure and 12,483 population controls from Italy, Spain, Norway and Germany/Austria, including stratified analyses based on age, sex and disease severity, as well as targeted analyses of chromosome Y haplotypes, the human leukocyte antigen (HLA) region and the SARS-CoV-2 peptidome. By inversion imputation, we traced a reported association at 17q21.31 to a highly pleiotropic ∼0.9-Mb inversion polymorphism and characterized the potential effects of the inversion in detail. Our data, together with the 5th release of summary statistics from the COVID-19 Host Genetics Initiative, also identified a new locus at 19q13.33, including NAPSA, a gene which is expressed primarily in alveolar cells responsible for gas exchange in the lung.Andre Franke and David Ellinghaus were supported by a grant from the German Federal Ministry of Education and Research (01KI20197), Andre Franke, David Ellinghaus and Frauke Degenhardt were supported by the Deutsche Forschungsgemeinschaft Cluster of Excellence “Precision Medicine in Chronic Inflammation” (EXC2167). David Ellinghaus was supported by the German Federal Ministry of Education and Research (BMBF) within the framework of the Computational Life Sciences funding concept (CompLS grant 031L0165). David Ellinghaus, Karina Banasik and Søren Brunak acknowledge the Novo Nordisk Foundation (grant NNF14CC0001 and NNF17OC0027594). Tobias L. Lenz, Ana Teles and Onur Özer were funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation), project numbers 279645989; 433116033; 437857095. Mareike Wendorff and Hesham ElAbd are supported by the German Research Foundation (DFG) through the Research Training Group 1743, "Genes, Environment and Inflammation". This project was supported by a Covid-19 grant from the German Federal Ministry of Education and Research (BMBF; ID: 01KI20197). Luca Valenti received funding from: Ricerca Finalizzata Ministero della Salute RF2016-02364358, Italian Ministry of Health ""CV PREVITAL – strategie di prevenzione primaria cardiovascolare primaria nella popolazione italiana; The European Union (EU) Programme Horizon 2020 (under grant agreement No. 777377) for the project LITMUS- and for the project ""REVEAL""; Fondazione IRCCS Ca' Granda ""Ricerca corrente"", Fondazione Sviluppo Ca' Granda ""Liver-BIBLE"" (PR-0391), Fondazione IRCCS Ca' Granda ""5permille"" ""COVID-19 Biobank"" (RC100017A). Andrea Biondi was supported by the grant from Fondazione Cariplo to Fondazione Tettamanti: "Biobanking of Covid-19 patient samples to support national and international research (Covid-Bank). This research was partly funded by a MIUR grant to the Department of Medical Sciences, under the program "Dipartimenti di Eccellenza 2018–2022". This study makes use of data generated by the GCAT-Genomes for Life. Cohort study of the Genomes of Catalonia, Fundació IGTP. IGTP is part of the CERCA Program / Generalitat de Catalunya. GCAT is supported by Acción de Dinamización del ISCIIIMINECO and the Ministry of Health of the Generalitat of Catalunya (ADE 10/00026); the Agència de Gestió d’Ajuts Universitaris i de Recerca (AGAUR) (2017-SGR 529). Marta Marquié received research funding from ant PI19/00335 Acción Estratégica en Salud, integrated in the Spanish National RDI Plan and financed by ISCIIISubdirección General de Evaluación and the Fondo Europeo de Desarrollo Regional (FEDER-Una manera de hacer Europa").Beatriz Cortes is supported by national grants PI18/01512. Xavier Farre is supported by VEIS project (001-P-001647) (cofunded by European Regional Development Fund (ERDF), “A way to build Europe”). Additional data included in this study was obtained in part by the COVICAT Study Group (Cohort Covid de Catalunya) supported by IsGlobal and IGTP, EIT COVID-19 Rapid Response activity 73A and SR20-01024 La Caixa Foundation. Antonio Julià and Sara Marsal were supported by the Spanish Ministry of Economy and Competitiveness (grant numbers: PSE-010000-2006-6 and IPT-010000-2010-36). Antonio Julià was also supported the by national grant PI17/00019 from the Acción Estratégica en Salud (ISCIII) and the FEDER. The Basque Biobank is a hospitalrelated platform that also involves all Osakidetza health centres, the Basque government's Department of Health and Onkologikoa, is operated by the Basque Foundation for Health Innovation and Research-BIOEF. Mario Cáceres received Grants BFU2016-77244-R and PID2019-107836RB-I00 funded by the Agencia Estatal de Investigación (AEI, Spain) and the European Regional Development Fund (FEDER, EU). Manuel Romero Gómez, Javier Ampuero Herrojo, Rocío Gallego Durán and Douglas Maya Miles are supported by the “Spanish Ministry of Economy, Innovation and Competition, the Instituto de Salud Carlos III” (PI19/01404, PI16/01842, PI19/00589, PI17/00535 and GLD19/00100), and by the Andalussian government (Proyectos Estratégicos-Fondos Feder PE-0451-2018, COVID-Premed, COVID GWAs). The position held by Itziar de Rojas Salarich is funded by grant FI20/00215, PFIS Contratos Predoctorales de Formación en Investigación en Salud. Enrique Calderón's team is supported by CIBER of Epidemiology and Public Health (CIBERESP), "Instituto de Salud Carlos III". Jan Cato Holter reports grants from Research Council of Norway grant no 312780 during the conduct of the study. Dr. Solligård: reports grants from Research Council of Norway grant no 312769. The BioMaterialBank Nord is supported by the German Center for Lung Research (DZL), Airway Research Center North (ARCN). The BioMaterialBank Nord is member of popgen 2.0 network (P2N). Philipp Koehler has received non-financial scientific grants from Miltenyi Biotec GmbH, Bergisch Gladbach, Germany, and the Cologne Excellence Cluster on Cellular Stress Responses in Aging-Associated Diseases, University of Cologne, Cologne, Germany. He is supported by the German Federal Ministry of Education and Research (BMBF).Oliver A. Cornely is supported by the German Federal Ministry of Research and Education and is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany's Excellence Strategy – CECAD, EXC 2030 – 390661388. The COMRI cohort is funded by Technical University of Munich, Munich, Germany. Genotyping was performed by the Genotyping laboratory of Institute for Molecular Medicine Finland FIMM Technology Centre, University of Helsinki. This work was supported by grants of the Rolf M. Schwiete Stiftung, the Saarland University, BMBF and The States of Saarland and Lower Saxony. Kerstin U. Ludwig is supported by the German Research Foundation (DFG, LU-1944/3-1). Genotyping for the BoSCO study is funded by the Institute of Human Genetics, University Hospital Bonn. Frank Hanses was supported by the Bavarian State Ministry for Science and Arts. Part of the genotyping was supported by a grant to Alfredo Ramirez from the German Federal Ministry of Education and Research (BMBF, grant: 01ED1619A, European Alzheimer DNA BioBank, EADB) within the context of the EU Joint Programme – Neurodegenerative Disease Research (JPND). Additional funding was derived from the German Research Foundation (DFG) grant: RA 1971/6-1 to Alfredo Ramirez. Philip Rosenstiel is supported by the DFG (CCGA Sequencing Centre and DFG ExC2167 PMI and by SH state funds for COVID19 research). Florian Tran is supported by the Clinician Scientist Program of the Deutsche Forschungsgemeinschaft Cluster of Excellence “Precision Medicine in Chronic Inflammation” (EXC2167). Christoph Lange and Jan Heyckendorf are supported by the German Center for Infection Research (DZIF). Thorsen Brenner, Marc M Berger, Oliver Witzke und Anke Hinney are supported by the Stiftung Universitätsmedizin Essen. Marialbert Acosta-Herrera was supported by Juan de la Cierva Incorporacion program, grant IJC2018-035131-I funded by MCIN/AEI/10.13039/501100011033. Eva C Schulte is supported by the Deutsche Forschungsgemeinschaft (DFG; SCHU 2419/2-1).N

High-speed exhaustive 3-locus interaction epistasis analysis on FPGAs

[Abstract]: Epistasis, the interaction between genes, has become a major topic in molecular and quantitative genetics. It is believed that these interactions play a significant role in genetic variations causing complex diseases. Several algorithms have been employed to detect pairwise interactions in genome-wide association studies (GWAS) but revealing higher order interactions remains a computationally challenging task. State of the art tools are not able to perform exhaustive search for all three-locus interactions in reasonable time even for relatively small input datasets. In this paper we present how a hardware-assisted design can solve this problem and provide fast, efficient and exhaustive third-order epistasis analysis with up-to-date FPGA technology.This study makes use of data generated by the Wellcome Trust Case–Control Consortium. A full list of the investigators who contributed to the generation of the data is available from http://www.wtccc.org.uk. Funding for the project was provided by the Wellcome Trust under award 076113 and 085475.Wellcome Trust; 076113Wellcome trust; 08547

Protein-coding variants contribute to the risk of atopic dermatitis and skin-specific gene expression.

Background: Fifteen percent of atopic dermatitis (AD) liability-scale heritability could be attributed to 31 susceptibility loci identified by using genome-wide association studies, with only 3 of them (IL13, IL-6 receptor [IL6R], and filaggrin [FLG]) resolved to protein-coding variants.Objective: We examined whether a significant portion of unexplained AD heritability is further explained by low-frequency and rare variants in the gene-coding sequence.Methods: We evaluated common, low-frequency, and rare protein-coding variants using exome chip and replication genotype data of 15,574 patients and 377,839 control subjects combined with whole-transcriptome data on lesional, nonlesional, and healthy skin samples of 27 patients and 38 control subjects.Results: An additional 12.56% (SE, 0.74%) of AD heritability is explained by rare protein-coding variation. We identified docking protein 2 (DOK2) and CD200 receptor 1 (CD200R1) as novel genome-wide significant susceptibility genes. Rare coding variants associated with AD are further enriched in 5 genes (IL-4 receptor [IL4R], IL13, Janus kinase 1 [JAK1], JAK2, and tyrosine kinase 2 [TYK2]) of the IL13 pathway, all of which are targets for novel systemic AD therapeutics. Multiomics-based network and RNA sequencing analysis revealed DOK2 as a central hub interacting with, among others, CD200R1, IL6R, and signal transducer and activator of transcription 3 (STAT3). Multitissue gene expression profile analysis for 53 tissue types from the Genotype-Tissue Expression project showed that disease-associated protein-coding variants exert their greatest effect in skin tissues.Conclusion: Our discoveries highlight a major role of rare coding variants in AD acting independently of common variants. Further extensive functional studies are required to detect all potential causal variants and to specify the contribution of the novel susceptibility genes DOK2 and CD200R1 to overall disease susceptibility

Protein-coding variants contribute to the risk of atopic dermatitis and skin-specific gene expression

BACKGROUND: 15% of atopic dermatitis liability-scale heritability could be attributed to 31 susceptibility loci identified by genome-wide association studies, with only three of them (IL13, IL6R, and FLG) resolved to protein-coding variants. OBJECTIVE: We examined whether a significant portion of unexplained atopic dermatitis heritability is further explained by low-frequency and rare variants in gene coding sequence. METHODS: We evaluated common, low-frequency and rare protein-coding variants using exome chip and replication genotype data of 15,574 patients and 377,839 controls, combined with whole transcriptome data on lesional, non-lesional and healthy skin samples of 27 patients and 38 controls. RESULTS: Additional 12.56% (s.e. 0.74%) of atopic dermatitis heritability is explained by rare protein-coding variation. We identified Docking protein 2 (DOK2) and CD200 Receptor 1 (CD200R1) as novel genome-wide significant susceptibility genes. Rare coding variants associated with atopic dermatitis are further enriched in five genes (IL4R, IL13, JAK1, JAK2, TYK2) of the IL13 pathway, all of which are targets for novel systemic atopic dermatitis therapeutics. Multiomics-based network and RNA-Sequencing analysis revealed DOK2 as a central hub interacting, among others, with CD200R1, IL6R and STAT3. Multi-tissue gene expression profile analysis for 53 tissue types from GTEx showed that disease-associated protein-coding variants exert their greatest effect in skin tissues. CONCLUSION: Our discoveries highlight a major role of rare coding variants in atopic dermatitis acting independently of common variants. Further extensive functional studies are required to detect all potential causal variants and to specify the contribution of novel susceptibility genes DOK2 and CD200R1 to overall disease susceptibility